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We study lattice SU(2) Yang-Mills theory with dimension d > 4. The model can be 
expressed as a (d— l)-dimensional 0(4) non-linear cr-model in a ci-dimensional heat bath. 
As is well known, the non-linear cr-model alone shows a phase transition. If the quark 
confinement is a consequence of absence of a phase transition for the Yang-Mills theory, 
then the fluctuations of the heat bath must destroy the long-range order of the non-linear 
ex-model. In order to clarify whether this is true, we replace the fluctuations of the heat 
bath with Gaussian random variables, and obtain a Langevin equation which yields the 
effective action of the non-linear cr-model through analyzing the Fokker-Planck equation. 
It turns out that the fluctuations indeed destroy the long-range order of the non-linear 
cr-model within a mean field approximation estimating a critical point, whereas for the 
corresponding U(l) gauge theory, the phase transition to the massless phase remains 
against the fluctuations. 
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1 Introduction 

We study Euclidean SU(2) Yang-Mills theory on the hypercubic lattice Z d with dimension 
d > 4. It is widely believed thato the gauge theory shows a quark confinement phase with 
a mass gap for all the values of the coupling in dimensions d = 4. On the other hand, the 
corresponding U(l) gauge theory in dimensions d — 4 is proven to show the existence of 
a deconfining transition to a massless phase [21 13] . Thus it is expected that there exists a 
crucial difference between SU(2) and U(l) gauge theories. 

In this paper, we explore the origin of this difference. For this purpose, we go back 
to the paper by Durhuus and Frohlich [4]. They showed that the cZ-dimensional Yang- 
Mills system can be interpreted as many (d — l)-dimensional non-linear cr-models which 
are stacked up in the d-th direction and coupled through (d — l)-dimensional external 
Yang-Mills fields! When we give our eye to one of the (d — l)-dimensional non-linear 
cr-models, the system can be interpreted as a (d — l)-dimensional non-linear cr-model in a 
c/-dimensional heat bath. When we turn off the interaction between the non-linear cr-model 
and the heat bath, the non-linear cr-model becomes the standard 0(4) non-linear cr-model 
because the gauge group SU(2) is homeomorphic to 3-sphere S 3 . As is well known, the 0(4) 
non-linear cr-model is proven to show a phase transition [7] in dimensions greater than or 
equal to three. This implies that, if the quark confinement is a consequence of absence of a 
phase transition for the Yang-Mills theory, then the fluctuations of the external Yang-Mills 
fields must destroy the long-range order of the 0(4) non-linear cr-model. 

The effective action of the (d — l)-dimensional non-linear cr-model can be derived 
by integrating out the degrees of freedom of the heat bath. However, carrying out the 
integration is very difficult. Instead of doing so, we replace the fluctuations of the external 
Yang-Mills fields with Gaussian random variables. Within this approximation, the spins 
of the non-linear cr-model can be interpreted as "particles" which move on S 3 , acted by 
the two-body interaction and the random forces. Namely the dynamics of the "particles" 
obeys a Langevin equation [5]. As is well known, a Langevin dynamics yields a Fokker- 
Planck equation which describes the time evolution of the distribution of the "particles" . 
In the present system, the effective action of the non-linear cr-model can be derived from 
the steady state solution to the corresponding Fokker-Planck equation. In the effective 
action so obtained, the attractive potential between the two "particles" is modified by the 
fluctuations of the external Yang-Mills fields. 

We show that the height and the width of the barrier of the attractive potential depend 
on the coupling constant of the Yang-Mills theory. Roughly speaking, the critical value 
of the coupling constant for the phase transition to a massless phase can be estimated by 
the height and the width of the barrier of the attractive potential. Therefore the critical 
value becomes a function of the coupling constant. In consequence, we obtain that within 
a certain mean field approximation, the critical value is always strictly less than the value 
of the coupling constant itself for weak couplings. This implies that the critical value must 
be equal to zero, i.e., there is no phase transition to a massless phase for non-zero coupling 
constants. 

On the other hand, the corresponding U(l) gauge theory shows that the attractive 



1 See, for example, the book pQ. 
2 See also related articles [SJH]. 
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potential does not depend on the coupling constant for weak coupling constants within 
the same approximation. Namely the fluctuations of the external Yang-Mills fields does 
not affect the critical behavior of the 0(2) non-linear cx-model. 

This paper is organized as follows. In the next section, we express SU(2) Yang-Mills 
theory in the form of the 0(4) non-linear er- model with a large heat bath, following 
Durhuus and Frohlich |4]. In Section [3j we obtain the Langevin equation for the "parti- 
cles" moving on § 3 , by replacing the fluctuations of the heat bath with Gaussian random 
variables. In the standard procedure, the Langevin equation yields the Fokker-Planck 
equation for the distribution of the "particles". In Section [U a steady state solution to the 
Fokker-Planck equation is obtained. The result immediately yields the effective action of 
the non- linear cx-model. Further, we show that the phase transition of the 0(4) non-linear 
cx-model disappears, owing to the fluctuations, within a mean field approximation for the 
effective action so obtained. In Section 0, we apply the same method to the correspond- 
ing U(l) gauge theory, and show that the phase transition to the massless phase remains 
against the fluctuations. 



2 Yang-Mills theory as a cr-model in a heat bath 

Let A be a sublattice of Z d . The SU(2) gauge field on A is a map from the oriented links or 
nearest neighbour pairs (q, q') of sites, q, q', of the lattice A into the Lie group G =SU(2), 

(q,q')^^eG, (2.1) 

obeying 

Oq'q = (tV)^ • (2-2) 

Let 7 be an oriented path which is written 7 = (q x , q2)(q2, qs) • • • (qn-ij Qn) with the 
oriented links, (q«,qi+i) of the neighboring sites, q^q^+i, for i = 1, 2, . . . , n — 1. When 
qi — q™, the path 7 is a loop. For an oriented path 7, we write 

U-y — ^qiq2^q2q;i ' ' * ^qn-iqn" (2-3) 

The Euclidean action of pure Yang-Mills theory on the lattice A C Z d is given by 

^ M (A):=-ij>eTrC/ ap , (2.4) 

where p denotes an oriented plaquette(unit square) of A, and dp is the oriented loop 
formed by the four sides of p. The orientation of the loop dp obeys the orientation of the 
plaquette p. The expectation value is given by 

(' " ->a := Z K l f J] dU b (- • •) exp [-/3^ M (A)] (2-5) 

J bcA 

with the inverse temperature /3 and the normalization where b is a link in A and dUf, 
is the Haar measure of the gauge group G =SU(2). 



Fluctuations in SU(2) Yang-Mills Theory 



4 



Following Durhuus and Frohlich [3], we use the relation between the d- dimensional 
Yang-Mills action and a (d — l)-dimensional non-linear a- model. The coordinates of a 
lattice site q are denoted (x^\x^ 2 \ . . . = (i,x^) with i = (x^\ . . . ,x^ d ~^) E 

7h d ~ x . Write A T = A fl {q : x^ = r} for the (d — l)-dimensional hyperplane, and A = 
A n 7h d ~ l x {0} for the projection onto 7h d ~ x lattice. Let (r) denote the gauge field £/ qq ' 
assigned to the link (q, q') in A T with q = (i, r) and q' = (j,t), and U[(r) the gauge 
field [/ qq ' with q = (i, r) and q' = (i, r + 1). The former are called horizontal gauge fields 
localized at 

x (d) = 

and the latter are called vertical gauge fields localized in the slice 
[r, r + 1]. Now the Yang-Mills action can be rewritten as 

^J M (A) = -^E ReTrf 4-^E E ReTrUyirr'U^Uyir^r+l). (2.6) 

T PC At T (ij)cA 

The first term in the right-hand side is a sum of Yang-Mills actions which depend on the 
horizontal gauge fields in (d — l)-dimensional hyperplane at x^ d ' = r. As to the second 
term, the vertical gauge fields in different slices are not coupled to each other. Therefore 
the summand about r in the second term is written in an action of a (d — l)-dimensional 
non-linear a-model for the vertical gauge fields as 

AU(^U h (r),U h (T + l)) = ~ J2 BeTtVr{T)- l U*{T)U?{T)u£(T + l) (2.7) 

(ij>CA° 

in the external gauge fields, U h (r) = {E/J(r)} and U h (r + 1) = {Ujfcr + 1)}. 

Let § 3 denote the 3-sphere. In order to express the gauge fields in terms of spins 
S G § 3 , we use the homeomorphism p : S 3 — > SU(2) which is defined by [4] 

„(S) = p (S™ 5« S«) = ( ^ + ) (2.8) 



with the radius {S^) 2 + {S^) 2 + {S^) 2 + {S^) 2 = 1. Then the interaction potential V l2 
between two spins Si and S2 in the non-linear a-model (12.71) can be written 

V 12 = -~ ReTr ^ (Si)" 1 p (a x ) p (S 2 ) p {a^ 1 , (2.9) 

where we have written <xi and <x 2 for the external horizontal gauge fields. When the 
external gauge fields, erg, take the vacuum configurations, cri = cr 2 = (1,0,0,0), the 
interaction becomes that of the 0(4) non-linear a- model in (d — 1) dimensions as 

1 3 
V 12 = -- ReTr p (Sx)" 1 p (S 2 ) = -Si • S 2 = - S[ k) S ( 2 k) . (2.10) 

k=0 

As is well known, the 0(4) non-linear a- model shows a long-range order of spins at low 
temperatures in three or higher dimensions [7] . The long-range order leads to the perimeter 
law of the decay of the Wilson loop [I] . The perimeter law implies deconfinement of quarks. 
If the confinement of quarks indeed occurs in the SU(2) gauge theory, the fluctuations of 
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the external gauge fields around the vacuum must destroy the long-range order of the O (4) 
non-linear cx-model. 

In order to take account of the fluctuations around the vacuum configuration of the 
external gauge fields, we approximate &£ as 



l-\&e\ 2 ,<Ti) ~(l,<x*) (2.11) 



with small fluctuations, 



r (i) J?) J$ 



for £ = 1,2. (2.12) 



We write 5cr£ = (0, Then the two-body potential is written 

V l2 « -Si • S 2 - 1 ReTr ^(S 1 )-V / (^i)^(S 2 ) - \ Re Tr ^(S^V^y (-5o- 2 ), (2.13) 
dropping the second ordei^l in the fluctuations b<Ji. Here we have written 



The right-hand side of (12.131) can be written 

V 12 ~ Vq + Vr (2.15) 



with 
and 
where 

and 



V = -Si • S 2 (2.16) 
Vr = -V2o-+ • (Si x S 2 ) - V2&- ■ (s{ 0) S 2 - 5f S x ) , (2.17) 

<x± = -^(<x 2 ±<xi), (2.18) 

8 e =(sJi L) ,S?\sl 3)> ), 1=1,2. (2.19) 

Thus the present system can be expressed as the 0(4) non-linear a-model in the heat 
bath. The interaction between the non-linear cr-model and the heat bath is given by Vr. 



3 The contributions of the second order of the fluctuations Scr? give order of temperature T = /3 _1 in 
the potential Vi2- Therefore one can expect that the contributions of the second order slightly modifies 
the coupling constants of the interaction potentials at low temperatures. 
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3 Langevin dynamics for two particles on S 3 . 

If we can integrate out the degrees of freedom of the heat bath, then we can obtain the 
effective action of the non-linear a-model. However, it is very difficult problem. Instead 
of this way, we replace the fluctuations of the external gauge fields with Gaussian random 
variables. Then, the spins of the a-model can be interpreted as the "particles" which move 
on § 3 , acted by the two-body interaction and the random forces. 

In order to derive the effective two-body interaction between two spins of the a-model 
within this approximation, we first introduce the Langevin equation for the two "particles" . 
We write xi = (xf' ,xf\xf'), £ = 1,2, for the local coordinates of the two 3-spheres § 3 . 
Then the Langevin equation [5] is given by 

j t xf = F$ + F^ £ = 1,2; . = 1,2,3. (3.1) 

with the forces, F^j, F^\, which are given by the gradientQ of the potentials as 

F§ = -g%*Vo (3.2) 

and 

= -9%dj,iV R , (3.3) 

where g % \ is the matrix inverse of the metric tensor g^^ for the "particle" £, and we have 
used the Einstein summation convention and written 

Let Pt{xi,x 2 ) be the distribution of the two "particles" on § 3 x § 3 . The expectation 
value of the function f{x\, x 2 ) on S 3 x S 3 at time t is given by 



(f)t '■= / f(xi,x 2 )pt(xi,x 2 )dfixdfi2, (3.5) 

where we have written 

dpi = y/ det gt dx^pdxf^dxf^ for £=1,2. (3.6) 
For a small At > 0, the following relation must hold: 

(f)t+At = ^ [ f{x l {t + At),x 2 {t + At))p t {x l ,x 2 )dp l dp 2 + 0{{Atf), (3.7) 

where E stands for the average over the fluctuations &e, £ = 1,2, and xe(t + At) is the 
solution of the Langevin equation ( 13. ip with the initial conditions xe(t) = xe at time t. As 



'See, for example, the book [5]. 
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usual, we assume that, for the short interval [t,t + At], the fluctuations frg are constant, 
and satisfy 



E 



a 







0, E 



a (i) a U) 



a 

AT 



5 lJ and E 



afaf 



a 
At 



S ij , 



(3. 



where a and a' are a nonnegative constant, and 8^ is the Kronecker delta. Physically, a 
natural assumption is that a and a' satisfy the condition a > a' > 0. From the relation 
between the fluctuations and the temperature of the heat bath, both of a and a' are 
proportional to the temperature of the heat bath. 
From the Langevin equation ( 13.11) . we have 



xf\t) 



i doc p i^ij 
dt 



dt 



dt'Ff\x{t')) 



(3.9) 



where we have written F^ 1 = F^j + F^' g and x(t) = (xi(t),X2{t)). Using this relation, we 
obtain 

jf>(x(0) = if } (^)) + E dF ^ t]) f dt " F ^")) + ■■■■ ( 3 - 10 ) 

™ t, OXm J t 



m,k 



Combining these, the expansion with respect to At is derived as 



x 



(0 



(* + At) = xfit) + F«(£(t))At + \ dF ^ t)] Fi k \mmy< 

OXm 



+ ■ ■ 



(3.11) 



m,k 



Substituting this into (13 .7p and using (I3.8p . the order of At yields 

d^ 



dpf{x)^— 



+ 



M 

At 

T 

At 



(JJb Q 



dfi 



d 2 f(x) 

-» /--b rv* ^ 



E 



m e ~ d dx\'dx) 



R,mA 



+ / ^ E 



0/(5;) 



(0 



E 



(3.12) 



where we have written M = § 3 x S 3 and d/x = dp\dp2- Since this equation holds for any 
function /, we can derive the equation of the time evolution for the distribution p t , i.e., 
the Fokker-Planck equation. 

To this end, consider first the first term in the right-hand side of (13.121) . Note that 



F$(£)p,(£) 







^— i ( k V det Ji F ( J (x) / (x)p t (x) 
Vdet & Q X W 



- £/(*) 



5 



Vdet g e F^j(x)p t 



x 



Vdet ^ a x « 

div £ [-F ,*(x)/(x)/>t(x)] - /(x) div £ [F ^(x)p t (x)] , (3.13) 
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where div^ stands for the divergence for the "particle" I. Combining this with the diver- 
gence theoremjf) 



dfx £ div e v e = 0, 



s 3 



(3.14) 



for a vector field Vi on S> 3 , the first term in the right-hand side of (13.1 2p is written as 

l d/i (duf) F®pt = ~J2 [ d vf div ^ ( F °JPt) ■ ( 3 - 15 ) 
ti Jm e Jm 

As to the second and third terms in the right-hand side of (13.121) , we must compute the 
second moments of the random forces. But one can treat these terms in the same way as 
in the above. The detail is given in Appendix |A] As a result, the Fokker-Planck equation 
is given by 



(a + a') {divi [r} x W ■ div 2 ('n 2 p t )} + div 2 [r] 2 W ■ divi^pt)]} 
2a'^div m C m ■ div„(C n Pt) 



(3.16) 



where is the Laplacian for the "particle" £, and we have written W = Si ■ S 2 ; the vector 
fields, £t, r] e and Ce, are given by 



C := g%*W, 



and 



7(0) « 



(3.17) 
(3.18) 

(3.19) 



for i = 1, 2, 3 and £ = 1,2. Here the vectors r\\ have four components like S^, and Ce have 
three components like S^. This Fokker-Planck equation can be written 



= — div J with div J = divi J\ + div 2 J 2 

at 

in terms of the current J = (Ji, J 2 ) which is given by 



P — n lJ T 



(3.20) 



(3.21) 



with 



idj,iV )pt - (a + a') {d jA p t - [(d jt l W)diy 1 ^ lPt ) + W(d jt A) ■ div 2 (r7 2 p t )]} 



-2a'Cj,i ■ divi(CiPt) + div 2 (C 2 p f ) 



(3.22) 



and with J j 2 given by exchanging the subscripts 1 and 2 in Jj X - Here we have written 



;) / q(0)c q(°)c 



(3.23) 



'See, for example, Theorem 5.11 in Chap. II of the book [5]. 
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4 A steady state for the Fokker-Planck dynamics 

The effective potential V e s between the two "particles" is derived from a steady distri- 
bution p t = p for the Fokker-Planck equation (13.201) . as in (14. 7p below. For a steady 
distribution p t = p, the Fokker-Planck equation (I3.20p becomes div J = 0. In order to 
obtain the solution near the north pole, = (1,0,0,0), for £ = 1,2, we introduce the 
local coordinates, (xg, ye, Zg) for £=1,2, as 



sjl-xj-yj- zj,x t , ye, z e ) . (4.1 



We write 

r = (x, y, z) = (xi - x 2 ,yi - y 2 , z x - z 2 ) (4.2) 

and 

R = (X, Y, Z) = (xi + x 2 , yi + y 2 , z x + z 2 ). (4.3) 

We also write r = |r| and R — |R|. In order to solve the partial differential equa- 
tion div J = 0, we employ the Cauchy-Kowalevski type expansion^! with respect to small 
x i, Uii z i- 

Let us compute the x-component J x> i of the current J\ for the particle 1. Note that 
V = -Si • S 2 = -1 + \r 2 + \{v ■ R) 2 + • • • . (4.4) 

/ o 

Immediately, 

^ = x + h r . R)x + h r . R)X + .... (4.5) 
ox\ 4 4 

Therefore, the first term of J Xt i of ( I3.22p becomes 



(d xA V )p 



-x - ~(r ■ K)x - ~(r ■ R)X + 



p. (4.6) 



In order to treat the rest of the terms of J Xt \, we assume that the steady state solution 
p t = p of div J = has the form, 

p = exp[-pV eS ], (4.7) 

where V e s is the effective potential to be determined, and (3 is the inverse temperature of 
the heat bath. Both of a and a 1 are proportional to the temperature as mentioned 
in the preceding section. The effective potential V e s must be vanishing for r = because 
the two-body potential (I2.13P becomes constant irrespective of the external fluctuations 
for Si = S 2 . From this and taking account of the spherical and exchange symmetries, we 
assume that the effective potential V e s can be expended as 

V eS = C 20 r 2 + C 40 r 4 + C 22 r 2 R 2 + C' 22 (r ■ R) 2 H , (4.8) 



'See, for example, Sec. D of Chap. 1 in the book [llj . 
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where C 20 , C 40 , C 22 and C 22 are the coefficients to be determined. In the following, we take 
a and a 1 to be small, and ignore the order of a and a'. 
For small xe, ye, zg, the current J x \ is written 



JxA 



-x - -(r ■ R)a; - -(r ■ R)X 



+ (a + a 



x x 



p — [a — a 
dp dp dp 



d d 



dx\ dx^ 



P 



dp dp dp \ r 2 dp 

+ y^— + z^- ) + x [ x 2 - — h y 2 ^ — h z 2 — 

dx2 dy 2 dz 2 



+ 2a' 



dx\ dyi dz\ 
dp dp dp 



2 dx? 



dx\ dy\ dz\ 



xi x- h y- h z—— + x 2 \ xi- h yi- h Zi^— 



dp 



dp dp 



dxi dyi 



dz\ 



dp 



dp dp 



- [ -x + -X ) x 2 - h y 2 - h z 2 — 



dxo 



'dy 2 



'dZ? 



, dp 

] dx~^2 yi1 



I dp 

dxo 



+ 



(4.9) 



The derivation is given in Appendix Q3] Let us substitute p of (14. 7p with the effective 
potential (14. 8 p into this right-hand side. First of all, since the leading order which is 
proportional to xexp[— f3V e s] must be vanishing, we have 



Since we can choose 



4/3 (a - a')C 20 = 1. 
1 







a — a' 



without loss of generality, we have 



a 



20 



Using these, we get 



with 



— (a — a 



d 



d 



d 



dxi dx 2 



exp [-^Kff] 



dV,* dV e 



cff 



dx± dx 2 



exp [-/3Kff] 



d 



dxi dx 2 
Moreover we have 



V cS = x + 8C 40 r 2 x + AC 22 R 2 x + 4C 22 (r • R)X + 



X-2 



(x 9 


d 


\ dxi 




d d 


d 


dx\ ^ l dyi 


dzi 


d d 


d 


dx 2 dy 2 


dz 2 



d 



2' 

lp r 2 _ . R) + 

l(3r 2 + ^(3(r-R) + 



exp[-f3V eS ] 
exp[-/3V eS ] 



r 2 dp 
~2~dx~ 2 



P 2a 
— — xr + ■ 

4 



exp[-/3Kff] 



4.10) 
4.11) 
4.12) 

4.13) 
4.14) 

4.15) 

4.16) 
4.17) 
4.18) 
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and 



- To 



dp 



+ o( r i + r 2 



dp (5 



dxi 2 K 1 z 'dx 2 4 
Substituting these into (14. 9D . we obtain 



[ r 2 + R 2 _ (r . R) ] exp[-/3Kff] + 



(4.19) 



J Xtl exp[/3V eS \ 



\8C, 



40 



+ 



4C 22 + 



1] r 2 x + 



R 2 x + 



[r 2 X - (r • R)x] 

a + a')j3 



4C 22 - 



From div J = 0, the coefficients must satisfy the relations, 

5(8^40 - 1) + ot'fi = 

and 



4C 2 2 + 



a'/3 



+ 



4C 22 - 



(a + a')(3 



0. 



; r -R)X + ---. (4.20) 

(4.21) 
(4.22) 



Using these relations, the current J x ^ can be written 



Jx,i = |-^r 2 a; + ^ [r 2 X - (r ■ R)x] + A [R 2 x - 3(r ■ R)X] \ exp[-f3V e{i ] + ■ ■ 



with the constant, 



A = 4C 22 + 



(4.23) 



(4.24) 



which we cannot determine in the present method. Clearly one notices that in div J, there 
appear the other terms, 



1 



a'fi 2 r A and - Af3[r R - 3(r • R 



(4.25) 



These are higher order in powers of the local coordinates but order of /3. Since the equation 
div J = must hold, this implies that there must exist some terms of order of (3 in the 
effective potential V e s so as to cancel the above terms of (14.251) . 

When both of the coefficients C 22 and C 22 depend on 0, the corresponding terms may 
appear in the expansion. In this case, from f)4.22p . we have 



C 22 ~ Cp and C 22 ~ -3C/3 
with some constant C for a large ft. Substituting these into V e s, we have 



~ -r 2 + C i0 r A + CP[r 2 R 2 — 3(r • R) 



(4.26) 



(4.27) 



This leads to instability of binding of the two particles because the value of R 2 is expected 
to become larger than order of in the thermal equilibrium. Thus we require that both 
of C 22 and C 22 are order of 1. 
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In consequence, we need the following terms in the effective potential V e ^\ 

C 60 r 6 , C i2 r A R\ C' 42 r 2 (r • R) 2 . (4.28) 

Here all the coefficients, C 60 , C 42 , C' 42l are proportional to (3 for a large (3. In the same way 
as in the above, we can determine these coefficients as 

C 6 o = -|a'/3 2 , C 42 = and C' A2 = (4.29) 

7! 5o 5o 

so as to cancel the above terms (I4.25P which appear in div J. As a result, the dominant 
contributions in the effective potential V e s are given by 

V cS ~ -r 2 - |a'/3 2 r 6 + ±-A{3r 2 \r 2 R 2 - 3(r • R) 2 ] (4.30) 
4 7! 5o 

for a large /3 because the second, third and fourth terms in the right-hand side of (I4.8P do 
not affect the critical behavior. 

Now we discuss the critical behavior of the (d— l)-dimensional a model with the above 
two-body interaction V e ^ . Consider first the case of A = 0. Namely the effective potential 
is given by 

V eS ~ \r 2 - |a'/3 2 r 6 (4.31) 

for small r and large 0. The second term lowers the potential barrier. Within a mean-field 
approximation [12], the critical temperature Tq can be estimated by the volume and the 
height of the potential well. More precisely, Tc ~ (volume) x (height). In the present case, 
the width w and the height h of the effective potential V e s are estimated as 

w ~ (A/3)- 1/4 , h ~ (X/3y 1/2 , (4.32) 

where we have written 

3' 

A = 12 • -a 13. (4.33) 

Therefore the critical temperature T c is estimated as 

T c ~ w 3 x h ~ (A/3)" 5/4 . (4.34) 

This is lower than for small temperature T = f3~ 1 . This implies that the true critical 
temperature must be equal to zero. 

In the case of A ^ 0, the third term in the right-hand side of (I4.30P may heighten 
the potential barrier if R 2 does not take a small value. But it is impossible that the 
term heightens the potential barrier in all the directions of r. Thus we reach the same 
conclusion, Tc = 0. 

Let us make the following two remarks: 

1. Our argument can be applied to the systems in arbitrary dimensions. Therefore a 
reader might think that our method suggests no phase transition for non-Abelian 
lattice gauge theory also in five or higher dimensions. On this point, we should 
remark the following: We used the two-body approximation, considering only a 
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single plaquette. When dealing with two plaquettes within our method, three- and 
four-body interactions would appear in the effective potential for the non-linear o- 
model. The resulting interactions may change the conclusion of this section. Namely 
a high- dimensional system may exhibit a phase transition. Actually, in five or higher 
dimensions, the effect of the three- or four-body interactions may not be ignored 
because the number of the neighboring plaquettes for a fixed plaquette becomes large, 
compared to low-dimensional systems. However, taking account of such interactions 
is not so easy. 

2. Consider the 0(4) non-linear a- model on the three-dimensional lattice with the 
effective two-body interaction which we obtained. Then the correlation length of 
the model leads to an estimate of the string tension [H [5] . Does the scaling limit so 
obtained give the standard continuum? This problem must be very important. But 
it is very difficult to compute the low temperature asymptotics of the correlation 
length for such a weakly attractive potential. 



5 Difference between U(l) and SU(2) gauge theories 

Let us see difference between U(l) and SU(2) gauge theories. 

For this purpose, we apply the present method to the abelian case G =U(1). In the 
case, the gauge field Ub on a link b is written 

U b = exp[i6 b ] (5.1) 

in terms of the angle variable B b G [0, 2tt). Therefore the two-body interaction V±2 between 
6>i and 62 is written 

V 12 = - cos(9 1 -02 + cn- <j 2 ), (5.2) 

where <j\ and <Ji are the angle variables of the external fields. We write 9 = Q\ — 9 2 and 
5a = <j\ — cr 2 , and assume that 5a is a small fluctuation. Under this assumption, the 
potential can be approximated as 

V12 ~ -cos6» + 5o-sin6>. (5.3) 

Then the Langevin equation is given by 



As usual, we assume 



d8 

— = — sin 9 — 5a cos 9. (5.4) 

dt y ' 

E[(H 2 ] = ^ (5.5) 

for a small At. In the same way as in the SU(2) case, we obtain the Fokker-Planck 
equation, 



dt 



2 



d ex d ex d 

-sin£ + --sin£cos£ + -^cos 2 £ 



Pt- (5.6) 
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For a steady state p t = p, we have 



• n a ■ n n a & 2 

sintH — sinflcoscH — tt cos i 
2 2d9 



p = 0. 



One can easily find the solution, 



(cos0) _1 exp [-2a- 1 /cos9], for -vr/2 < 9 < vr/2; 



0. 



otherwise. 



(5.7) 



(5i 



Since the diffusion disappears at 9 = ±7r/2 in the right-hand side of (j5.4p . the "particle" 
cannot move beyond the points. Clearly, we have 



p ~ const. exp[— a l 9 2 ^ 



(5.9) 



for a small Thus there is no term which is proportional to a^ 1 or higher powers of a~ l 
in the effective potential, and the critical behavior can be expected to be the same as the 
standard 0(2) nonlinear-cr model. This is consistent with the rigorous result of [21 [3]. 



A Derivation of the Fokker-Planck equation 

Consider first the case with a' = in (13.81) . We introduce o l i satisfying = — cr 1 - 7 with 



(JOl 02 03\ 

[a ,<j ,<j J 



and (<7 23 ,cr 31 ,<7 12 ; 
Then the random potential Vr of (12.171) can be written 

\r 1 _ „ij o(k) q W 

VR — 7= £ ijk£ CT Oi 2 , 

V2 



(A.l) 



(A.2) 



where e^t is completely antisymmetric and satisfies £0123 — +1, and we have used the 
Einstein summation convention. From a' = 0, we have 



E [a af3 a mn ] 



a 
At 



(A.3) 



Using (1A.2|) and (1A.3j) . we obtain 
E[(d £tl Vn) (d kjl V R )\ 



2 

a 



2At 

2a 



At 



£ [(dtA" ] ) (d k A J) ) S { 2 5] S? - (d u S[ 7) ) 4 7) (d k A 5) ) S?] • (A.4) 

7,5 



Using the metric 



9ij,e 



dx^ dx 



(A.5) 
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of § 3 for the "particle" £, the above result is written 



and 



In 

E [(d u V R ) (0 M V H )] = ^ \gtk,i - (d e ,iW)(d kjl W)\ 

2a 

E [(9^Vk) (d k , 2 V R )] = — [g ek ,2 - (d e , 2 W)(d k!2 W)} , 



where we have written W = S\ ■ S 2 . Similarly, we have 



E[(4,AiVr) (dt,iV K )] 



Combining this with 



2a 
At 



7,<5 



h) 



q (S) q (S) 
J 2 J 2 



fl2n(l) QCiW 

c ^1 0(7) 
dx^dx^p dxf^ 



f)2q(l) a oil) 



we obtain 



2a 



E[(d kA d jA V R ) (d iA V R )] = — - (^.ifyiWO (Qj.iW)] 



where 1 are the Christoffel symbols [9]. In the same way, we get 



and 



E [(d e>1 V R ) (d k ,2V R )} = ~W (d e>1 d k>2 W) 

2a 

E \{d K2 d hl V R ) (d ifl V K )] = (9 fc)2 W) (fy !^W0 • 



(A.6) 
(A.7) 

(A.8) 
(A.9) 

(A.10) 

(A.ll) 
(A.12) 



Using ( 1A.6I) . we have 
E 



R.l R.l 



E 
2a 



g\ (dt,iV R ) cf k (d k>1 v R ) 
g u i9 j i[9i k ,i-{d t ,xW){d ktl w)] 



At 



(A.13) 



where Q is the vector field which is given by f)3.17p . From (1A. 6|) and (lA.lOj) . we obtain 



k r n,i 



+ 



E [{d k , ig \d hl V R ) {g k {d^V R )] 

(d k J{) g k [E [(d jA V R ) (d eA V R )} + g\g^i E [(d ki idj t iV R ) (d £A V R )} 
( Wi) 9 k i \9it,i - {di,iW) (d^W)} 

2a 

At 
2a 



— [d hl g\ + g'iTl^ - {d k ^[) g] 



(A.14) 
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where we have used3 



(A.15) 



In the same way, the relations ( 1A.11I) and ( 1A.12I) yield 



E 



ei(0 p(J) 

R,l R,2 



and 



cte* 



E [(W^iVk) (s^Vr)] 



^E[(^, 2 %i^r) (dtjVn)) 



(A.16) 



(A.17) 



respectively. The contribution from the two random forces Fr^ with the same indexes 
I = 1 in the right-hand side of ( I3.12p becomes 



'11 



E 



^i^'ll + £ / ^ 



» A: 



5/ 



E 



dx^ 



F (k) 

R.l 



Pt, 



(A.18) 



where we have used (1A. 13|) and (1A.14j) . Note that 

d 2 f 



M 



Pt 



diif(Axpt) 



(A.19) 



where the second equality follows from the property^! of the Laplacian A^. The rest of the 
contributions in the right-hand side of (1A.18I) are computed as 



M 



df 



Pt 



7 See, for example, Sec. 7 of Chap. I of the book [TO] . 

8 See, for example, Corollary 5.13 in Chap. II of the book [H). 
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r ["J i 

= / d V Tx^ ^ 1 V^9i( d 3,if)CAlfh - (d jt if) == ^,i\/det^^ijPf 

= - I dp 7=^= — fyi v 7 det fli&Pt 
= - d/j (d jA f)£i div! [£ x p f ] 

</Af 

= - f dfi 1 <9j,i-v/detffi£j/div x (giPi) 
7 M V det 5-i 

+ / rift / ~== 9j,i V det i^j divi (^ipt) 
y M Vdet 5-i 

d/i/divi^idiviCCiPt)], (A.20) 



'AT 

where we have used the divergence theorem (13.141) . Substituting this and (1 A. 19f) into 
(jA.18jl . we obtain 

In = af ^/{A lPt -di Vl [6di Vl (eip t )]}. (A.21) 

Next consider the contribution from the two random forces with different indexes, 
£ = 1 and £ = 2, in the right-hand side of (13. 12ft . Using ( 1A.16I) and ( 1A.17I) . we obtain 



'12 



/ * E TTJT9'i9"(*.='' | -')(^,i*,2lV)p 1 



P (fc) 



-a 



-a 



dfi J— djjy/det 9~2g j 2(di,if)gi(di,id k:2 W)W Pt 
m Vdet 5(2 



+a [ d^(d i> J)^=d jt2 ^^ 2 g j ^g u 1 (d e , 1 d ki2 W)W Pt 
J M Vdet g 2 

-a [ dfi(d h J)g\g k ^d k , 2 W)(d^d e , 2 W)p t 

JM 

= a! ^(^x^^W^^^V^et^^'^^i^^W)^, (A.22) 
7 M Vdet 5f 2 

where we have used the divergence theorem (I3.14p . Recalling W = Si • S2, we have 

9/,i^,2W = (^,iSi) • (<9 fc , 2 S 2 ) . (A.23) 
Substituting this into the above result, we get 

/ dp {d i , 1 f)g i {(d t , 1 S 1 )W ■ -7^=d ji2 y/tet^tfi (<9 fc , 2 S 2 ) p t 
Jm Vdet g 2 



'12 = ot 
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a / dfj, {d it if)rf x W ■ div 2 {^Pt) 



M 



a \ dfj, = — di^&etgxifJW ■ &iv 2 {r)2pt) 
Jm V det #1 

-a / d/x/ = ^,iV^et^>yiW • div 2 {r} 2 Pt) 

J M Vdets-i 

-a / <i/i / divi [77^ • div 2 (r) 2 Pt)] , 
Jm 



(A.24) 



where rf t is given by (13181) . From fl3~T2l . (13~T5]) . (|AJ8|) . ( TA~2Tj) . f[A22]) and ([Oil] , we 
obtain the Fokker-Planck equation, 



~ir = div£ ( F WP^ + a ^2 i^ePt - div £ div^(&p t )]} 

-a {div! [tjjW • div 2 (?7 2 pi)] + div 2 • div^ry^)]} 



(A.25) 



for a' = 0. 

Next consider the case with a' ^ 0. To begin with, we note that 



E 



(0 0') 



E 



- ^ E[a«a?] + E[^Vi J; ] + E[<W] + Eft% 



Oh + <7i 



«„Cj)i 



a + a' 



5 ij . 



Similarly, 

Further, we have 
E 



E 



a — a 



At 



(A.26) 



(A.27) 



(*) 0') 



-E 



0. 



Since we can write 



E 



(*) U) 
a a 



a + a ' 5 n _ ^L 8 io 

At At ' 



(A.28) 



(A.29) 



it is sufficient to calculate the corrections from the second term in this right-hand side, 
with replacing a with a + a' in the above result (1A.25I) . 

In (1A.13I) . the correction to E ^(c^iVr,)^ (c^iVr) is given by 



4a' ~i p 



(A.30) 
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where Ce 1S given by (I3.19p . Similarly, the correction to E [(dk ! ig l {dj i iVYi)(g k (de : iV^)~\ in 
flA.141) is given by 

4a' / „ »j \ ik , . * 

(A.31) 



Therefore the same calculations as those from flA. 18|) to (1A.21j) yield the correction, 



— 2a'divi 



Ci ■ div i(CiP^ 



(A.32) 



in the right-hand side of the Fokker-Planck equation (1A.25|) . 

In f ]A.16j) . the correction to E g l x{d^{VjC)g^ 2 {%, , 2Yk) is given by 

_ 4a' p ~j 
At 

Further, the correction to E [(dk^idj^V^^g^de^VR)] in flA. 17[) is given by 

- % (dU\) ■ & 

Therefore similar calculations to those from ( 1A.22I) to ( 1A.24I) yield the correction, 



(A.33) 



- 2a'divi 



(A.34) 



(A.35) 



in the right-hand side of the Fokker-Planck equation (1A.25I) . In consequence, the Fokker- 
Planck equation is given by 



dpt 
dt 



- ^ div ^ ( F o,ePt) + (a + a') ^ {^iPt ~ div ^ [& ^i(^ept)]} 

t i 
-(a + a') {div! {q^W ■ div 2 (r7 2 p t )] + div 2 [q 2 W ■ divifopt)]} 

-2a' dW m Cm ■ dWn(CnPt) 



(A.36) 



B Derivation of the expansion (14.91) 

The metric g^j of S 3 is computed as 

/ 1 + j£xj jmye lixizt \ / 1 + xj x e yi 

gij,e= ItV&i t + leVe im*i = Ve x e 1 + ViH | H , (B.l) 

\ nztxt Itzm ^ + lezjj \ zex e z e y £ 

where we have written 




H = , = with n = Jxj + yf + z^. (B.2) 
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Therefore, the inverse g\ is given by 



9 i 



Using this, we have 

(0x,iWOdivi(£ip) 



l-xj -x e y e -x e z e 
-y e x e 1 - y\ -y e z £ | - 
-Z(X t -z e y e 1-2$ 



-^^(^iSi • S 2 )d iA p H 

dp dp dp \ 

x 1 x x — *~y* — ^ Z TT~ i + 

ox i oyi oz\ 



Similarly, 

W{d x ^x) ■ div 2 (r7 2 p) = W{d x ^x) • ^(^ 2 S 2 )9 i)2 p + 

= ^(^■ 2 s 1 -s 2 yi^ 2 p+-- 

H ^,2 -x-^(r.R)x! + . 



0%2P + 



1 



<9p <9p <9p 

2^27; h ^i?/2q h ^1^2^— 

0x2 dy 2 dz 2 



+ 



1 - -r' 



„nd^ 2idp_ 31V 
y 2 o 2 ^y 2 o 
cte 2 ay 2 dz 2 . 



+ 



dp dp dp 

dx 2 dy 2 dz 2 



+ 



dp _ dp 



dx 2 2 ' dx 2 
dp 



dp 



2 dp dp 
dx 2 dy 2 dz 2 



+ 
dp 



dp dp 

xix 2 —— + x ± y 2 - h xiz 2 —— 

dx 2 dy 2 dz 2 



+ 



1 2 dp 
— —r — h x 



We write 
Note that 



dx 2 2 dx 2 

> _ /7(1) A(2) a(3)\ 



(9/9 (9/9 <9p 

hy 2 ^ ^2^— 

dx 2 dy 2 dz 2 



+ 



Ad) _ 

Sx,l ~~ 



si a) - 



1 — r 



2 an 



Therefore, we have 

Cx,l 



-X1X2 



1 - r 2 2 , 



-X1Z2 



(B.3) 



(B.4) 



(B.5) 
(B.6) 



(B.7) 
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-xix 2 - \ 1 - 72, -Xi?/2, -a;i^2 



In the same way, 



-2/i^2, -ViVz 



1 - r|, -yi22 + • 



and 



C,i = \ -ZiX 2 , -ziy 2 , -ziz 2 - yjl-rl \ - 
From these results, we obtain 

C,i • L,i = l-rj + 2x x x 2 + • • 

Cx,i • C„,i = z/i^a + + ■■■ 

CxA ' CzA = z \ x 2 + x l z 2 H 



and 



(B.8) 

(B.9) 
(B.10) 

(B.H) 
(B.12) 

(B.13) 



Using these, we have 
Cr.i -divi(CiP) 



= Cx,i • 9 l \Cj,A,ip + ■■■ 

= (1 - r 2 2 + 2x l x 2 )g i \d i>1 p 

+ (yi^ 2 + x\y 2 )g % \di±p + (>ix 2 + x^g^d^p + 



(l-ri + 2x 1 x 2 ) 



dp 



dp 

xiyi- xxzi 



2 



dp dp 



— r 2 — x\ I x— h y— h 2 



<9yi 
f- ■ 

dp 



dp 
dz\ 



dp dp dp 

+ x 2 [xi- \-yx- h zi — 

dx\ dyi ozi 



dx\ dyi dz\ 



In the same way, 



C x ,2 = \ xix 2 + y 1 - rf , x 2 yi, x 2 ^i J H . 

C y ,2 = (y2Xi,yiy 2 + yjl-rj, y 2 zij +■■■ 



and 



C 2 , 2 = ( z 2 xi,z 2 y 1: z\z 2 + y 1 - r{ ] + 
Combining these, (1B.8|) . (1B.9j) and ( IB.lOj) . we obtain 



Cx,i ■ Cx,2 = - ( 1 - \r\ ~ \r\ + 2xix 2 ^ + 



(B.14) 

(B.15) 
(B.16) 

(B.17) 
(B.18) 
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and 



L,i-C y ,2 = -2xm + --- (B.19) 
Li ■ L2 = -2*1*2 + • • • • (B.20) 



Using these, we have 

Cx,i • div 2(C 2 p) = Cx,i • fotjlPiflP + 



- - \r\ ~ \r\ + 2*1X2) g^d^p 

-IxxyigV^izP - 2x 1 z 2 g l \d i ^p H 

-g^p+^rl + rl)^- 

( dp dp dp \ 

-2xi x 2 - V Vi q V z 2 ^— + ■ 

\ ox 2 oy 2 oz 2 J 



dx 2 2 dx 2 



Substituting (jMD, f lRij) . flR5|) . flRlij) and f lR2T|) into (ET22]L we obtain the expansion 

dMD. 
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